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Particle distributions in squeezed states, even and odd coherent states are 
• given in terms of multivariable Hermite polynomials. The Q-function and 

\0 • Wigner function for nonclassical field states are discussed. 
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^ : 1 Introduction 

The aim of the talk is to give a review of such nonclassical states of light as squeezed 
^ . states IH, 0, correlated states even and odd coherent states (Schròdinger 

! cat states 0), squeezed Schròdinger cat states 0. First we discuss the photon dis- 

qh| tribution function for the generalized correlated states of multimode light. For 

finding Solutions of Schròdinger equation with time-dependent Hamiltonians which 
■ are genèric quadratic forms in position and momentum operators the integrals of mo- 

'jjj . tion which are linear forms have been constructed in Refs. ^ and |TD|. Such 

I integrals of motion have been analyzed and applied to general problems of quantum 

mechanics and statistics in Refs. ||11||,||12||. 



The nonstationary Hamiltonians are appropriate models for the physical condi- 
tions in which nonclassical states of fields (photons, phonons, pions, etc.) may be 
created. Following ||T3| we will give the result for photon distribution function of 
genèric mixed squeezed and correlated Gaussian state. Initially the state is taken 
to be Standard coherent one (partial case of such state is photon vacuum) and, for 
example, due to nonstationary Casimir effect it becomes multimode mixed correlated 
state. 



2 Multimode Hermite Polynomials and Mixed Cor- 
related Light 

The most general mixed squeezed state of the N-mode light with a Gaussian density 
operator q is described by the Wigner function W^(p, q) of the genèric Gaussian 
form, 

^(Q- < Q >)M-^(Q- < Q > 



VF(p,q) = (detM)"2 exp 



1) 



where 2N-dimensional vector Q = (p, q) consists of N components pi, . . . , and 
N components qi, çat, operators p and q being the quadrature components of 

the photon creation àf and annihilation à operators (we use dimensionless variables 
and assume ^ = 1): 



P 

q 



àt 
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2N parameters < Pi > and < Qí >, i = 1,2, 
are the average vàlues of the quadratures, 



(2) 



A^, combined into vector < Q>, 



< P > 

< q > 



Tr ^p, 
Tr ^q. 
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A real symmetric dispersion matrix M consists of 2N^+N variances 
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2 \QaQf3 + QisQc 
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1,2,. ..,2N. 
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They obey certain constraints, which are nothing but the generalized uncertainty 
relations [|T2| . 



The photon distribution function of this state 



Pn = Tr^|n><n|, n = (ni, «2, • • • , «Ar), (5) 
where the state |n > is photon number state, which was calculated in [Q, and it is 



n! 



(6) 



The function Hj^^{y) is multidimensional Hermite polynomial. The probability to 
have no photons is 
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det (m + ilajv) 


exp 



< Q > (2M + hNT^ < Q > 



(7) 



where we introduced the matrix 

R = 2U^(1 + 2M)-^V* - aNx, (8) 

and the matrix 

In 



The argument of Hermite poljTiomial is 

y = 2U*(l2jv-2M)-i<Q>, (10) 
and the 2N-dimensional unitary matrix 

u = ^f-f^ (11) 

is introduced, in which I^r is the NxN identity matrix. Also we use the notation 

n! = ni\n2\...ni\f\. 

The mean photon number for j-th mode is expressed in terms of photon quadrature 
means and dispersions 

< >= \{^Pm + ^m. - 1) + >' + < >)• (12) 

We introduce a complex 2N~vector B = /?2, ■ ■ ■ ^ I3n, l3l, 13^, . . . , jd*^). Then 
the Q-function fï^ is the diagonal matrix element of the density operator in the 



coherent state basis | /32, . . . , /Sat > . This function is the generating function 
for the matrix elements of the density operator in the Fock basis |n> which has 
been calculated in [T^. In notations corresponding to the Wigner function (1) the 
Q-function is 

Q(B) = Po exp[-^B(i? + aNx)^ + Bi?y]. (13) 

Thus, if the Wigner function (1) is given one has the Q-function. Also, if one has the 
Q-function (13), i. e. the matrix R and the vector y, the Wigner function may be 
obtained due to the relations 

M = U*(i? + a^,)"'U^-^, 
<Q> = U*[l-(i? + a^,)-Vjv.]y. (14) 

For pure squeezed and correlated state with the wave function 

^ = A^exp[— xmx + cx], (15) 



where 

N = [det(m + m*)]^/^-^/^ exp{^(c + c*)(m + m*)-\c + c*)}, (16) 

the symmetric 2Nx2N-matrix R determining Q-function has the block-diagonal 
form 

(IT) 

The NxN-matrix r is expressed in terms of symmetric matrix m 

r* = l-{m + l/2)~\ (18) 
and the 2N-vector y = (Y, Y*) is given by the relation 

Y* = i=(m-l/2)-V (19) 

The corresponding blocks of the dispersion matrix 

M=( ^P'í ] (20) 

are 



(jqq = Í(m + m*)-\ 

o"pq = -{m - m*){m + m*)'-^. (21) 

The probabihty to have no photons is 

^ _ [det(m + m*)]^/2 
° ~ |det(m + 1/2)1 

« exp{l(c + c·)i,n + ™r '(c + c.) H- i|c(,„ + l/2)-c + c'i,n' + l/2)-'c1}. (22) 

The multivariable Hermite polynomials describe the photon distribution function for 
the multimode mixed and pure correlated hght |TB|, I^.The nonclassical state 



of the hght may be created due to nonstationary Casimir effect |JÏ8|], and the Husimi 
oscillator is the model to describe the behaviour of the squeezed and correlated pho- 
tons. 



3 Multimode Even and Odd Coherent States 



We define the multimode even and odd coherent states (Schròdinger cat male states 
and Schròdinger cat female states, respectively) as 



I A± >= iV±(| A > ± I -A >), (23) 

where the multimode coherent state | A > is 

I A >=| «1, «2, . . . , a„ >= -D(A) I >, (24) 

and the multimode coherent state is created from multimode vacuum state | > 
by the multimode displacement operator -D(A). The definition of multimode even 
and odd coherent states is the obvious gener alization of the one-mode even and odd 



coherent state given in |1TT|]. The normalization constants for multimode even and 
odd coherent states become 

|A|2 

e 2 

= 



2./cosh I A P' 



N- = / (25) 
2A/sinh A P 



where A = ( «i, «2, • • • , «n) is a complex vector and its modulus is 

n 

I A P = | «1 |2 + I «2 r + • . . + I «n r= I] I 1^ • (26) 



m=l 



Such multimode even and odd coherent states can be decomposed into multimode 
number states as 



A± >= iV± V , ' 'U^" (1 ± (-l)-i+«2+--") I n >, (27) 



where | n >= | rzi, n2, . . . , n„ > is the multimode number state. Also from Eq. 
(23), we can derive an important relation for the multimode even and odd coherent 
states, namely. 



a,- I Al > = an/tanh I A p I A_ >, 



ai I A_ > = an/coth | A P | A+ > . 

(28) 



The probability of finding n photons in multimode even and odd coherent states can 
be worked out with the help of Eq. (27) 



«1 


2ni 


«2 


2n2 


1 '-'íji 


|2n„ 


(ni!)(n2! 
1 «1 


)... 

«2 


(n„!) cosh . 

1 o- 1 


2n„ 



[nil)[n2í) . . . (n„!) sinh | A p 



(29) 



Multimode coherent states are the product of independent coherent states of each 
mode, and photon distribution function is the product of independent Poissonian 
distribution functions. But in the present case of multimode even and odd coherent 
states we cannot factorize their multimode photon distribution functions due to the 
presence of the nonfactorizable cosh | A p and sinh | A p. This fact implies the 
phenomenon of statistical dependences of different modes of these states on each 
other. 

In order to describe the properties of the distribution functions from Eq. (29) we 
will calculate the symmetric 2Nx2N dispersion matrix for multimode field quadrature 
components. For even and odd coherent states we have 

< A± I aiük I A± >= «jOífc, (30) 

and complex conjugate vàlues of the above equation for < A± \ ajal \ A± >. Since 
the quantity < A± | | A± > is equal to zero the above equation gives two NxN 
blocks of the dispersion matrix. For other two NxN blocks of this matrix we have 
for the multimode even coherent states 

1 1 
(7+t =< A+ I -(atofc + akül) I A+ >= a-Ofetanh | A ^ +-6ik, (31) 

and for multimode odd coherent states 

(7~t , A_ I i(ajafc + akül) \ A_ >= a*ak coth | A +^Sik. (32) 
z z 

For the dispersion matrix, the mean vàlues of the photon numbers rii — a|aj for 
multimode even and odd coherent states are the foUowing 



nj A+ 
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p tanh 


1 A 1 
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> = 


ai ^ coth 


A|' 



< Ah 

< A_ I I A_ > = \ ai\'^ coth | A |^ . (33) 



Taking into account the above equation the symmetric NxN dispersion matrices for 
photon number operators can be obtained from the above given distribution functions 
for multimode even and odd coherent states. By defining 



afk =< A± I mnk I, A± >, (34) 



the corresponding expressions in such states are 



on 
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P sech ^ 



A|^ + 



ik 



«i I I «fc 



csch ^ I A 



ai 1^ tanh | A 
f I «i P coth I A p Sik- 



(35) 



As the nondiagonal matrix elements of the dispersion density matrix are not equal 
to zero so we can predict that different modes of these states are correlated with 
each other. In other words, as we have mentioned before, there exist some statistical 
dependences of different modes on each other. 

Another interesting property for the multimode even and odd coherent states is 
the Q-function, and it can be obtained in the foUowing manner. First of all the 
density matrices for the multimode even and odd coherent states are 



A+ >< Aj 



(36) 



then the Q-function can be calculated as 



g+(B,B*) 

Q_(B,B*) 



< B 



< B I p_ 



B > 

! + |B|2) 

B > 

! + |B|2) 



cosh(AB*) 



sinh(AB* 



(37) 



where | B >=| /32, . . . , /3„ > is another multimode coherent state with multimode 
eigenvalue B = (/5i, /32, . . . , /^n)- We call these functions for even and odd coherent 
states as the Q-functions for the Schròdinger cat states. The Q-function for single- 
mode odd coherent state shows the cràter type behaviour for small vàlues of the 
quantity | a \ and for its larger vàlues the Q-function begins to split into two peaks 
in a similar manner as in case of even coherent states ÍIT^l. 



The Wigner function for the multimode coherent states is |]T2 
2^ exp[-2ZZ* + 2AZ* + 2B*Z - AB* 



W. 



A 



B 



A,B 



where 



q + ip 

V2 ■ 



(38) 
(39) 



For even and odd coherent states the Wigner function is 

W^A±(q, p) = \N±\^ [l^(A,B=A)(q, p) ± W^(A,B=-A)(q, p) 
± Vr(_A,B=A)(q, P) + W^(-A,B=-A)(q, p)], 



(40) 



where the explícit forms of N± are given in Eq. (25). For multimode case we use the 
foUowing notat ions 



AZ* = «iZ* + «2^2 + • ■ • "n^^, 

ZZ* = ZiZl + Z2Z; + . . . + ZnZ* 



(41) 



The photon distribution function gives the probability of finding 2k photons for 
two-mode even coherent state, and is defined as 
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(42) 



where 2k = ni + «2, for both rii and n2 to be even or odd numbers. Similarly for 
two-mode odd coherent state it gives the probabihty of finding 2k+l photons 



P_(2A; + 1) = 



(1 «1 




«2 


2-j2fc' 


fl 




(2fc + l)!sinh(| 


«1 


2 + 


«2 





(43) 



For this case, rii is even and «2 is odd number, or vice versa. For single-mode case 
the photon distribution functions demonstrate súper and sub-Poissonian properties 
for even and odd coherent states, respectively. The same conclusions may be drawn 
for two-mode (and multimode) even and odd coherent states. 

4 Some Relations for Wigner and Q— functions 

The Wigner function of a system W^(p, q) = W^(Q) is expressed in terms of density 
matrix in coordinate representation as (see, for example, |T2[) 

W{p, q) = y p(q + ^, q - ^) exp(-zpu) du. (44) 

The inverse transform is 

P(x, x') = J p) exp[zp(x - x')] dp. (45) 

The Q-function is expressed in terms of the Wigner function through the 3N-dimensional 
integral transform 

Q(B) = I <í>B(x,x',p)Py(^^,p) rfx í/x' dp (46) 

with the kernel 

$b(x, x', p) = 7r-^/2 exp[ïp(x - x') - ^B(or^, + l27v)B - Ç + v^Ba^v.X], (47) 

where the 2N-vector X = (x, x') is introduced. If one has the Q-function the Wigner 
function is given by the integral transform 

1 r ^ ~ ~ ~ 

W^(p,q) = ^ {\{d''(3k Afc du,^k{uk^B)]Q{B), (48) 
^ •' k=i 



where the argument of the Q-function B is replaced by the 2N-vector with complex 
components 

B = ( /32, . . . , (3n, 7Í, 72) • • • > 7iv)> 
and the kernel has the form 

$,K,B) = vr-V2exp[-|/3,p-|7,P + v^(g, + |^)/3, + v^fe-f)7| 



\{^^ + y)^ - ^(^. - y)^ - f - 1^ - W + 7,/3*]. (49) 



The density matrix in coordinate representation is related to the Q-function 

p(x,x') = n-''' J{Y[d'P, rfS0.(B) exp[--(4 + x,')]}QiB), (50) 
fc=i 

where the kernel of the transform is 

M^) = vr"'/' exp[-|/5fc|2 - |7fc|2 + v^x,/3, + V2xiY, _ ^ _ 2f + ^,^5*]. (51) 

The evolution of the Wigner function and Q-function for systems with quadratic 
Hamiltonians for any state is given by the following prescription. Given the Wigner 
function Vr(p, q, í = 0) for the initial time í = 0. Then the Wigner function for the 
time t is obtained by the replacement 

W{p,q,t)= W{p{t), q(í), í = 0), 

where the time-dependent arguments are the linear integrals of motion of the quadratic 
system found in |T2|. The same ansatz is used for the Q-function. Namely, given 



the Q-function of the quadratic system Q(B, {t = 0)) for the initial time í = 0. 
Then the Q-function for the time t is given by the replacement 

Q(B, t)= g(B(í), t = 0), 

where the 2N- vector B(í) is the integral of motion linear in annihilation and creation 
operators found in [jT2|. This ansatz foUows from the statement that the density 
operator of the Hamiltonian system is the integral of motion, and its matrix elements 
in any basis must depend on the appropriate integrals of motion. In particular, the 
Wigner function and Q-function depend just on the linear invariants found in 



5 Multivariable Hermite Polynomials 

For parametric forced oscillator the transition amplitude between its energy levels has 
been calculated as overlap integral of two genèric Hermite polynomials with a Gaus- 
sian function (Frank-Condon factor) and expressed in terms of Hermite polynomials 



of two variables ||^. For N-mode parametric oscillator the analogous amplitude has 
been expressed in terms of Hermite polynomials of 2N variables, i. e. the overlap 
integral of two genèric Hermite polynomials of N variables with a Gaussian function 
(Frank-Condon factor for a polyatomic molecule) has been evaluated in The 
corresponding result uses the formula 

{n = ni,n2, . . .riN, m = mi,m2, ■ . . rriN, m^, rij = 0, 1, . . .), 

/ iíí^>(x)iíM(Ax + d)exp(-xmx + cx) = ^==exp{-cm-'c)H^l{y), (52) 

V det m 4 

where the symmetric 2Nx2N-matrix 

with NxN-blocks Ri, R2, Ru is expressed in terms of symmetric NxN-matrices 
R, r, m and NxN-matrix A in the form 

Ri = R- -Rm-^R, 
2 

R2 = r — -rAm~^Ar, 

R12 = ~rAm-^R. (54) 

Here the matrix A is transposed matrix A and R12 is transposed matrix R12. The 
2N-vector y is expressed in terms of N-vectors c and d in the form 

(55, 

where the N-vectors yi and y2 are 

yi = ^{Rm~^ + m~^R)c 

y2 = -{rAm~^ + m'^Ar)c + rd. (56) 
For matrices R = 2, r = 2 the above formula (52) yields 

/ {Y[ HnXxi)H„^XYl ^ikXk + di)} exp(-xmx + cx) rfx 

i=l k=l 



-== exp(-cm-ic)iJÍ,^2(y), (57) 
vdetm 4 



with NxN-blocks Ri, R2, Ru expressed in terms of N x N-matrices m and A in 
the form 

R2 = 2{1- Am-^À), 

R12 = -2Am-\ (58) 
The 2N-vector y is expressed in terms of N-vectors c and d in the form (55) with 

yi = m"^c, 

ya = ^(Am"^ +m-^À)c + 2d. (59) 
If the symmetric matrix p has the block-diagonal structure 

P = 



Ri 

i?2 



with the symmetric SxS-matrix Ri and the symmetric (2N-S)x (2N-S)-matrix R2 
the multivariable Hermite polynomial is represented as the product of two Hermite 
polynomials depending on S and 2N-S variables, respectively, 

Ht\y)=HÍfHy,)HÍ^l{y2), 

where the 2N-vector y has the vector-components 

y = ( yi, y2), 

and 2N-vector k has the components 

k = (ns, n2Ar_s) = { ni, ns, ns+i, n2N)- 

The partial case of this relation is the relation for the Hermite polynomials with the 
matrix R with complex conjugate blocks i?i = r, R2 = r*, and complex conjugate 
vector-components yi = y2 

HfHy) = \HÍf\y,)\\ S = N. 

The calculated integrals are important to evaluate the Green function or density ma- 
trix for the systems with quadratic Hamiltonians. The partial cases of multivariable 



Hermite polynomials determine some other special functions [|T6|, |117[ . 



6 Squeezing in Parametric Oscillator 

For the parametric oscillator with the Hamiltonian 

where we take h = m = ci;(0) = 1, there exists the time-dependent integral of 
motion 

A = ^[e{t)p-é{t)x], (61) 

where 

e{t) + uj^it)e{t) = 0, £(0) = 1, 5(0) = i, (62) 
satisfying the commutation relation 

[A, A]] = 1. (63) 

It is easy to show that the packet solutions of the Schròdinger equation may be 
introduced and interpreted as coherent states since they are eigenstates of the 
operator A (61), of the form 

^,(x,t) = *o(a;,t)exp{-^^ + ^— }, (64) 

*o(x,t) = 7r-V4,(í)-i/2,^p!Í^ (65) 

is analog of the ground state of the oscillator and a is a complex number. The 
variances of the position and momentum of the parametric oscillator in the state (65) 
are 

2 ' ^ 2 

and the correlation coefficient of the position and momentum has the value corre- 
sponding to minimization of the Schròdinger uncertainty relation pi 



where 



= — - — , 0-p = — - — , (66) 



Another normalized solution to the Schròdinger equation 

^o.m{x, t) = 2iV^^o(a:, t) exp{-^^ cosh (68) 

where 

iv^^^ïHiaa (69) 

2A/cosh lap 



is the even coherent state (the Schròdinger cat male state). The odd coherent 
state of the parametric oscillator (Schròdinger cat female state) 



<í>^f{x,t)=2Nf<ÍIoix,t)exp{-L·J A^jsmh^-^, (70) 

N, ^ ^-^^tm, (71) 

2A/sinh lap 



where 



satisfies the Schròdinger equation and is the eigenstate of the integral of motion 
(as well as the even coherent state) with the eigenvalue a^. These states are one- 
mode examples of squeezed and correlated Schròdinger cat states constructed in 0. 



7 Conclusion 

The discussed nonclassical states of the fields (photons, phonons, pions) may be 
created in nonhnear interactions. The particle statistics with squeezing and corre- 
lations of quadrature components may give an experimental evidence of producing 
the new types of the field states. The Schròdinger-like equations are used also in 
other branches of physics like fiber òptics. In [^] the Schròdinger-like equations 
has been introduced to describe the charged particle beam in accelarator. The ap- 
proach described in the talk may be also applied in the classical physics using the 
quantum-mechanical methods. 



References 

[1] J.N. Hollenhorst, Phys. Rev. D19 (1979) 1669. 
[2] H.P. Yuen, Phys. Rev. A13 (1976) 2226. 

[3] V.V. Dodonov, E.V. Kurmyshev, and V.I. Man'ko, Phys. Lett. A79 (1980) 150. 
[4] V.V. Dodonov, LA. Malkin, and V.I. Man'ko, Physica 72 (1974) 597. 
[5] B. Yurke and D. Stoler, Phys. Rev. Lett. 57 (1986) 13. 

[6] V. V. Dodonov, V. I. Man'ko, and D. E. Nikonov, Even and odd coherent states 
(Schròdinger cat states) for multimode parametric systems, Preprint INFN-IY- 
NA-93/49, DSF-T-93/49 (1993); Phys. Rev. A (in press). 

[7] E. C. G. Sudarshan, Uncertainty relations, zero point energy and the linear 
canonical group, Proc. of Second International Workshop on Squeezed States 
and Uncertainty Relations, Moscow, 25-29 May, 1992, eds. D. Han, Y. S. Kim, 
and V.I.Man'ko (NASA Conference Pubhcation 3219, 1993) p. 241. 



[8] 1. A. Malkin and V. 1. Man'ko, Phys. Lett. A32 (1970) 243. 

[9] I. A. Malkin, V. I. Man'ko, and D. A. Trifonov, Phys. Lett. A30 (1969) 414. 

[10] I. A. Malkin, V. I. Man'ko, and D. A. Trifonov, Phys. Rev. D2 (1970) 1371. 

[11] I. A. Malkin and V. 1. Man'ko, Dynamical Symmetries and Coherent States of 
Quantum Systems (Nauka Publishcrs, Moscow, 1979) [in Russian]. 

[12] V. V. Dodonov and V. I. Man'ko, Invariants and Evolution of Nonstationary 
Quantum Systems, Proc. of Lebedev Physics Institute 183, ed. M. A. Markov 
(Nova Science, Commack, N. Y., 1989). 

[13] V. V. Dodonov, O. V. Man'ko, and V. I. Man'ko, Phys. Rev. A50 (1994) 813. 

[14] K. Husimi, Proc. Phys. Math. Soc. Japan 22 (1940) 264. 

[15] V. V. Dodonov, V. I. Man'ko, and V. V. Semjonov, Nuovo Cim. B83 (1984) 145. 
[16] V. V. Dodonov and V. I. Man'ko, J. Math. Phys. 35 (1994) 4277. 
[17] V. V. Dodonov, J. Phys. A27 (1994) 6191. 

[18] V. I. Man'ko, J. of Soviet Làser Research (Plenum Press, N. Y.) 12 (1991) N5. 

[19] N. A. Ansari and V. I. Man'ko, Phys. Rev. A50 (1994) 1942. 

[20] I. A. Malkin, V. I. Man'ko, and D. A. Trifonov, J. Math. Phys. 14 (1973) 576. 

[21] E. Schròdinger, Ber. Kgl. Akad. Wiss. Berlin 24 (1930) 296. 

[22] R. Fedele and G. Miele, Nuovo Cim. D13 (1991) 1527. 



